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On the Role of Network Centrality in the 
Controllability of Complex Networks 

Nicoletta Bof, Giacomo Baggio, and Sandro Zampieri 


Abstract —In recent years complex networks have gained in¬ 
creasing attention in different fields of science and engineering. 
The problem of controlling these networks is an interesting and 
challenging problem to investigate. In this paper we look at the 
controllability problem focusing on the energy needed for the 
control. Precisely not only we want to analyze whether a network 
can be controlled, but we also want to establish whether the 
control can be performed using a limited amount of energy. 
We restrict our study to Irreducible and (marginally) stable 
networks and we find that the leading right and left eigenvectors 
of the network matrix play a crucial role in this analysis. 
Interestingly, our results suggest the existence of a connection 
between controllability and network centrality, a well-known 
concept in network science. In case the network is reversible, the 
latter connection involves the PageRank, an extensively studied 
type of centrality measure. Finally, the proposed results are 
applied to examples concerning random graphs. 

Index Terms —Complex networks, controllability, network cen¬ 
trality, PageRank. 


1. Introduction 

Complex networks are systems composed of a large number 
of units which interact among themselves, forming in this 
way a behavior which is much richer than the behavior of 
the single units [1]. Many systems, which model both natural 
processes and engineering structures, can be seen as complex 
networks. Among them, one can mention genomic networks 
and ecologic networks in biology, social networks in sociology 
and economic or financial networks in economics, while in 
engineering, electric power grids, transportation networks and 
communication networks are some important examples [2]. 
Therefore many areas of science and technology can take a 
great advantage from a deep understanding of this class of 
systems. 

Consequently, the properties of complex networks have at¬ 
tracted a lot of interest among different scientific communities 
in the last years. Controllability is one of these properties, and 
it consists in the possibility of steering the state of the network 
from any initial value into a final arbitrary one, by fixing 
the profile of the state of a subset of nodes, called control 
nodes, which are assumed to be directly accessible by the 
controller [3]. This property has been analyzed by relating 
it to the features of the underlying graph [4]. In particular, 
the results on structural controllability theory [5] have been 
exploited. This kind of approach aims to establish whether the 
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network is controllable and, specifically, which nodes must be 
controlled in order to obtain the network controllability. 

From [4], several other papers focusing on this type of prob¬ 
lem have been proposed, including [6], [7], [8], [9]. However, 
an important aspect is neglected in this line of research. In 
fact, one can notice that, even if a control profile able to drive 
the network state may exist in principle, this profile may not 
be physically implementable due to the energy it requires. This 
observation gave rise to another type of approach related to 
network controllability, focusing specifically on the evaluation 
of a “degree” of controllability of the network. This concept 
can be made precise only by introducing a suitable metric 
related to the energy needed for the control. The evaluation 
of the control energy naturally involves, for linear systems, 
the notion of controllability Gramian, which is a symmetric 
positive definite matrix in case the network is controllable [3]. 
The analysis of the controllability Gramian offers different 
ways to study the degree of controllability of a network, 
depending on which property of this matrix is inspected. 
For instance, one can select the minimum eigenvalue of the 
controllability Gramian [10], [11], the trace of its inverse [12], 
its determinant [13], or its condition number [14]. 

In this context, different types of problems can be addressed. 
For example, one interesting issue consists in the optimal 
placement of the control nodes in a way such that the 
degree of controllability is maximized [12], [15], [16], [17]. 
Another engaging problem is to relate the network structure 
to the number of control nodes needed to make the network 
practically controllable, namely controllable by control profiles 
with bounded energy [10], [14], [11], [18], [19], [20]. 

In this paper, we focus on this second problem. Precisely, 
we study in which cases the energy needed for the control 
tends to infinity as the number of nodes gets larger and 
larger. This can be a way to establish the classes of complex 
networks which are practically impossible to control, since 
the energy they require for their control can go over any 
threshold as the number of nodes increases. This behavior 
will depend both on the properties of the network and on 
the number of control nodes. As in [11] and [18], we adopt 
the minimum eigenvalue of the controllability Gramian as the 
measure of the controllability degree. In our analysis we will 
restrict to networks with non-negative weight matrices which 
are irreducible and (marginally) stable. 

This type of analysis has been started in [18], where partial 
results have been proposed showing that isotropic networks are 
difficult to be controlled, while anisotropic networks are more 
controllable. However, only intuitive definitions of isotropic 
and anisotropic networks were given there, supported by 
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some illustrative examples. Actually, isotropic networks were 
described as networks in which there are no global preferential 
direction resulting from the network weights. Elaborating this 
idea, in the present paper we show that the right and left 
eigenvectors associated with the largest eigenvalue of the 
network adjacency matrix play an important role to determine 
whether the network is difficult to control. The importance of 
these eigenvectors is already well-known in network science, 
in connection with the concept of centrality, see [2, Section 
7.2]. Precisely, we show that when these eigenvectors suggest 
that all the nodes have similar centrality degrees, then the 
network will be difficult to control. 

Paper structure. The present paper is organized as follows. 
Section II is devoted to present some preliminary results. In 
Section III we state the main result of the paper. Section V 
contains some examples of application of the main theorem. 
Finally, in Section VI we summarize the contributions of the 
paper and we outline a number of possible future research 
directions. Most of the proofs of the more technical results 
are postponed to the appendix. 

Notation. We denote by M" the set of n dimensional vectors 
with real entries. The symbols and 1„ stand for the fc-th 
vector of the canonical basis and the n-dimensional vector with 
all ones, respectively. Moreover, we denote by the set 

of n X m matrices with real entries. The symbol /„ stands for 
the n X n identity matrix. Given a matrix A. the symbol Aij 
means the {i,j)-th entry of A, while A^ and ker(A) mean 
the transpose and the kernel of A, respectively. Finally, the 
symbol ||A|j 2 denotes the 2-norm of A. 

A matrix A € is said to be (Schur) stable if |jA ‘||2 —> 

0 as f —> oo, while it is said to be (Schur) marginally stable if 
11 A* II 2 is bounded in t. We denote the i-th component of the 
column vector v G ffi" as Vi. Given two vectors v,w G K", 
with u _L w, we indicate that v is orthogonal to w in the 
standard Euclidean metric, that is v^w = 0. We let diag(u), 
V G M", denote the n x n diagonal matrix with elements 
vi,... ,Vn on the diagonal. Given a positive definite matrix 
W G we indicate with \\x\\w '■= ^x^Wx , x G M", 

the weighted Euclidean norm (in case VE = /„ we simply 
write ||a:||). 

A matrix A is said to be non-negative (positive) if Aij > 0 
{Aij > 0) for all i,j = l,...,n. Irreducible and primi¬ 
tive matrices are special subclasses of non-negative matrices. 
Specifically, a non-negative matrix A G is said to be 

irreducible if {In + A)’^ is positive for some k G N, and 
primitive if is positive for some k gN [21, Chapter 8]. 

Also, we denote hy Q = (V, E, A) the weighted graph with 
vertex (or node) set )2 = {1,2,..., n}, edge set £ C V x V and 
adjacency matrix A G satisfying A^ > 0 if (j, i) G £. 

For a weighted graph Q = {V,£,A) the out-degree of the 
node i is equal to while the in-degree of the node 

i is equal to 'Y^jAij. A graph where the adjacency matrix 
has only 0-1 elements will be called unweighted. If A = A^, 
the graph is called undirected. In case of undirected graphs, 
the in-degrees are equal to the out-degrees and they are called 
simply degrees and denoted by the symbol deg^(i). In case 
the graph is undirected and unweighted, deg^(i) corresponds 
to the number of edges of the node i. 


Finally, given two functions f{x) and g{x), with g{x) non¬ 
zero, we write f{x) = o{g{x)) to mean that 


lim 

rc—>-oo 


f{x) 

9{x) 


= 0 . 


Other standard notation is taken from [11]. 


II. Preliminary facts 

Given a non-negative matrix A G the following 

assumption will be used throughout the paper. 

Assumption 1: Matrix A is irreducible and marginally stable. 

These conditions on A assure that there exists a unique real 
eigenvalue Ai such that 1 > Ai > |Ai| for all i = 2,..., n, 
where Ai denotes an eigenvalue of A; moreover, there exist 
positive left and right eigenvectors related to Ai [21, Theorem 
8.4.4]. 

One class of matrices which satisfies Assumption 1 are 
irreducible column-stochastic (row-stochastic) matrices, i.e. 
irreducible matrices whose columns (rows) sum up to 1. This 
is an important class which has a large number of applications, 
e.g. in the context of distributed optimization and consensus 
algorithms [22], [23], [24]. 

We denote with v and w the right and left eigenvector of 
A with respect to Ai respectively, that is 

Av = Aiv, A^ w = Aiw. (1) 


The vectors v and w are called the right and left leading 
eigenvectors of the matrix A. Since v and w are positive, we 
can define the following positive vector 


Wi 

W2 

Vl 

V2 


w. 
Vn J 


-,T 


( 2 ) 


and also the diagonal matrix 


n := diag(7r). 


( 3 ) 


The vectors v and w can be normalized in such a way that 
= 1- 

We exploit matrix If to define an analogous to the time 
reversal of an irreducible Markov chain, see e.g. [25, Section 
1.6], that is we introduce the following matrix 

A^ := U-^A^U. (4) 


This matrix is irreducible and v and w are respectively its right 
and left eigenvectors with respect to Ai. We can now define 
the matrix 

A^ .= n^/^AA^n-i/^ 

which is symmetric and positive semidefinite (and therefore 
has all real non-negative eigenvalues). Notice that, even if A 
and A^ are irreducible, AA^ (and therefore also A'^) may not 
be irreducible. Fet 


0-1 > 0-2 > • • • > (T„ > 0, (5) 

be the eigenvalues of A^'. The following Femma, whose proof 
is postponed to Appendix A, gives a characterization of the 
eigenvalues of A^. 
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Lemma 1: The largest eigenvalue of is equal to the 
square of the largest eigenvalue of A, namely ai = Af.' 

The next Proposition, whose proof is given in Appendix A, 
is instrumental for Section III. 

Proposition 1: If y G M" is such that y^v = 0, then for 
any f > 0 it holds 

l!(^^)‘2/lln-i<^2ll2/lln-i- 


Finding upper bounds for Aniin(VV’K:,T) can be therefore an 
instrument for proving that a system, although controllable 
in theory, is not controllable in practice due to the energy 
requirements. The following Theorem provides a useful bound. 

Theorem 1: Let A be a non-negative matrix which satisfies 
Assumption 1 and let the non-negative constant CT 2 and the 
positive vectors v and w be defined as in (5) and in (1), 
respectively. If the product AA^ is primitive, then 


III. The main result 

Suppose that we are given a network represented by the 
graph Q = (V,£,A), with A G satisfying Assumption 

1^. We consider the following discrete state-space system built 
according to the adjacency matrix of the graph; 

x(t + 1) = Ax{t) + Bfcutcit), 

where x G M", Bjc ■= • • • Ck^] G and ujc G M™, 

with 1C := {ki,... ,km} C V being the set of control nodes 
and m := \JC\. 

Our aim is to analyze the controllability of the system. 
However, we are not totally satisfied with the standard notion 
of controllability (from the origin) [3]: 

Definition 1 (Controllability): A system is controllable if, 
for every state Xf G M", there exists an input sequence {u(f)} 
that steers the initial state a;(0) = 0 to x/. 

In fact, we would also like to determine whether a control¬ 
lable system is difficult to control, where our idea of difficulty 
concerns the amount of energy needed to steer the initial state 
to a desired final one. To achieve this goal we introduce the 
controllability Gramian of the system, which is defined as 

T-l 

Wk.t ■.= Y,A^BkBI{A^)P 

T —0 

It is well known [3] that the system is controllable if and 
only if Wjc^T is invertible for a big enough T. For controllable 
systems, denote by m^(-) the unique minimum -energy 
input which steers the state from the initial value x(0) = 0 to 
the final value x{T) = Xf. The energy of u^(-) is given by 
(see [3]) 

T-l 

E{xf,T) := ^ ||w^(r)f = 

r—0 

From this observation it is possible to select a controllability 
metric. Indeed, choosing a worst-case analysis, we observe that 
the final state requiring the maximum energy is the one parallel 
to the eigenvector of 'Wjc,t corresponding to its minimum 
eigenvalue, which is denoted by Amin(VV’K:,T)- Precisely, 

max E(a;/,r) = Amin(W/c.T)“^ 
lk/ll=i 

In this way, we can say that a controllable system is dif¬ 
ficult to control if Aniin(VV’K:,T)~^ is big, or equivalently if 
Amin(W/c,T) is Small. 

^Indeed this result holds for any irreducible A, independently from its 
stability. 

^ We recall that irreducibility of A implies that the graph Q is connected. 


Amin(H’K:,T) < 


maxi{wi/vi} a 2 


(6) 


m\ni{wi/vi) (t|(1 - (T 2 ) 
for all T G N and all sets JC with cardinality m. 

Proof: Given T < T, the following chain of inequalities 
holds: 

Amin(W;(; r) = Hjin 


x\\ = l 


T-l 


= min x^ 'Wk,t - V' (A^)'^ x 

11 - 11=1 V / 

/ T-l 

> Amin(WK:.T) + min {-'Y' Bk:B^(A^)" 

lkll=i \ —- 


t=T 
/T-l 


= Amin(WK:,T) “ max x^ 

k =1 


Y^a^b^bKa^ 


\t^T 


We obtain in this way 

-^min (W^c.t) < A min (Wk.t) 


/T-l 


max X 
11-11=1 


YA^B^bI(A^)^ X. 


t=T 


Let us now define 


Amin(W’K:.T) := min x Wk,tx, 

11-11=1 

xJ-V 

where v is as in equation (1). As before, if T < T, we can 
prove that 


■^min (W^c.t) < A min (Wk,t) 


.T 




+ max X 

xXv \r—T / 

and moreover it holds that Amin(W’K:,T) < Amin(H’K:,T)- Tak¬ 
ing T = we have that Amin(yV’;c.T) = 0- Indeed, W^^.T 

coincides with where Cf jq is the controllability 

matrix of the system, i.e. 

Ct,k-~\.Bic ABjc A^B/c ■■■ A^~^Bic]- 

Now, if we introduce the matrix := [Cj, ^ u] then the 
rank of Ct^jcCj^ is strictly less than n. This in turn implies 
that there exists a x G M" such that x^CtxCjic^ = 0 and 
hence x is orthogonal to v and is such that x^W/c^tX = 0. 
From the fact that Aniin(H’;c t) = 0, we have 

-^min (Wk.t) < A min (Wk.t) 

< max x^ (A'^B!cB^(A^)'^'\ x < fif O) 


11-11=1 

xJ-V 


\t=T 
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where 

Bf := lim max 
T-)-oo 11 x 11 = 1 

XJ-V 

Our aim is to find a bound for /?j.. First observe that if we 
define a := min^ {vi/wi], due to the particular form of Bjq, 
it holds that 

y^BicBj-y < y^y < -y^R-^y = 

a a 

Using Proposition 1 and letting y = {A^yx, we obtain 

x^A-B^BZiA^Yx < 

a a 

In this way we obtain 
1 

Bt < — max ||x||n-i lim > cr?. 

a ||x|| = l T^oo ^ 

Av r^T 

Now, since ||x|| = 1, we have that 


A^BkBZ{A'Y X 


r=T 


ll^ll 


2 

n-i 


E Vi 

HI: 

2=1 


Vi 2 

xf 

Wj 


< 


max 
2=1,...,n 



max 
2=1,...,n 



Since A is marginally stable and AA^ is primitive, also A^ 
is primitive (since it has the same sparsity pattern of AA^). 
Therefore, it holds that a 2 is strictly less than ai = Xf < 1, 
and therefore is strictly less than 1. Summing up all previous 
considerations, we have 


uiaxYvi/wi} T maxiiwi/vi} crj 

—^ = lim > CT 2 = —^—f--. 

T[Vmi{VilWi\ T->oo Ti\\ni{Wi/Vi\ 1 — (72 

t—T 


Letting T = [^J, inequality (7) yields Xynini^KA < Pt- 
Since > n/m — 2 and < 1, the statement follows. 


Remark 1: Concerning the primitivity condition on AA^, 
a simple sufficient condition ensuring this property is A 
irreducible and with all its diagonal elements strictly greater 
than zero. As a matter of fact, if Aij > 0, then 

fAA^)ij = AiiAji + • • • + AijAjj + • • • + AinAjn > 0 . 

>0 

This entails that all the elements greater than 0 in A are 
greater than 0 also in AA^, therefore AA^ is irreducible and 
with positive diagonal elements, and thus primitive. This non¬ 
zero diagonal hypothesis is usually met by the networks we 
consider, since it implies that every node has its own dynamics. 

Remark 2: In case A is stable (and irreducible), the previous 
theorem holds tme even without requiring the primitivity of 
AA^. As a matter of fact, stability implies that Ai < 1 and it 
automatically follows that (J 2 < cri < 1. It is worth noting that 
this version of our result is stronger than what it is proven in 
[11, Theorem 3.1], namely that if 11A11 2 < 1, then the system is 
certainly difficult to control. As a matter of fact a stable matrix 
does not necessarily have 2-norm less than one. For instance. 


for all e such that 0 < £ < the matrix Ag := is 

stable, but ||Ae ||2 > 1. 

Once the system matrix A and the number of controlled 
nodes m are fixed, a metric that describes the controllability 
degree of the system is 


A(A,to):= inax Xmin{yVK,T)- 

/C : I /C I =Tn 
TeN 

Comparing several systems, those with a smaller A(A, m) are 
more difficult to control, that is they need more energy for 
their control. The previous theorem shows that 


A(A, m) < 


u\ayii{wi/vi) 
Tami{wi/vi} ct|(1 - ( 72 ) 


(8) 


Now, consider a sequence of systems with matrices A„ 
of increasing dimension n and a number of controlled nodes 
m(n) depending on n. In case 


A(A„,TO(n)) 0, (9) 


we have that, even if for all n the systems are controllable, for 
large n the energy required to control them is so high that the 
control is practically impossible^. When the right-hand side 
of (8) tends to 0 as n goes to infinity, (9) is verified and so 
the systems become practically uncontrollable when n is big. 
Therefore, if we want to have a chance that the systems are 
practically controllable, we need to fix m{n) in a way that the 
right-hand side of (6) does not go to 0 as n increases (even 
though, since we have only an upper bound, this fact alone 
does not guarantee that the practically controllability can be 
in fact achieved). 

To avoid that the right-hand side of (8) tends to 0 as n goes 
to infinity, we can choose a fast enough growing m{n). The 
growth of m{n) will depend on two characteristics of A: 


1) On how CT 2 depends on n: in some cases, as we will 
see in the next section, CT 2 stays bounded away from 
1, namely there exists a constant e > 0 independent 
of n such that 1 — (72 > e, for all n. In other cases, 
instead, (72 tends to 1 as n goes to infinity. In this case 
it is important to understand how fast this occurs. The 
quantity 1 — ( 72 , called the spectral gap of A, is essential 
to study inequality (8) in this case. 

2) On how the fraction _ max,{7r.} ^jgpgj^^jg 

mini {lUi/Ui I minilTTij ^ 

on n: as we will see in the following, in many important 
cases this fraction can be evaluated and in some of them 
this quantity remains bounded in n. 


The way these two sources of dependence on n interact will 
determine the growth of m{n). For instance, if (72 is bounded 
away from 1 and j-gmains bounded in n, then 

practical controllability can be achieved only if the number 
of control nodes m{n) grows linearly in n, namely only if a 
fixed fraction of the nodes are controlled. For example, this 
holds tme for symmetric matrices [18]. In fact, in this case, 
mtY{Y/Yj ~ following Proposition shows that this 

property holds for a more general class of matrices. 


^Note that the tightness of bound (6) is not relevant when we consider the 
asymptotic behaviour as n tends to infinity. 
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Proposition 2: Let v and w be the right and left eigenvector 
of A, then 




Proof: We have the following chain of equivalences 


max,{wi/uj 

T[mii{wi/vi} 


1 <J4> V = aw, a € M 

A^v = Aiv 
^ Av = Al^V 
^ {A- A^)v = 0 

<=> V € ker(A — A^). 


The class of matrices satisfying (10) comprises, for instance, 
irreducible doubly stochastic matrices, i.e. matrices which are 
both column- and row-stochastic. 


IV. Stochastic matrices. Eigenvector Centrality 
AND Reversible Matrices 


• If the stochastic matrix A is reversible, then A can be 
regarded as the transition matrix of an irreducible time- 
reversible Markov chain as defined in [25, Section 1.6]. 

• If 2l is reversible then it can be symmetrized by a diagonal 
transformation. Indeed, from (4) it follows that 

s := 

is symmetric, where 11 is defined in (3). This in turn im¬ 
plies that all the eigenvalues of A are real, and therefore, 
by virtue of the irreducibility of A, they can be ordered 
in decreasing order as follows 

1 = Ai > A 2 > • • • > A„ > —Ai. 

• For a reversible matrix A, we have that As = S'^ and so, 
besides having cti = Af as stated in Lemma 1, we also 
have that 

(72 = max{A2, A^}. (12) 

It is possible to associate with a given non-negative matrix 
C a stochastic matrix which is reversible if C is symmetric. 

Precisely, given an irreducible matrix C G we can 

obtain a column-stochastic irreducible matrix A as 


The result of Theorem 1 acquires a nice interpretation when 
matrix A is column-stochastic. Applying the theorem to such 
matrices, the next Corollary immediately follows: 

Corollary 1: If A is an irreducible column-stochastic matrix 
such that AAA is primitive, then for all T € N and all sets JC 
with cardinality m it holds that 


Amin(H’K:,T) < 


^max ^2 
Vmin 0 'i(l - < 72 ) 


( 11 ) 


with Umax := maxi Vi and u^in := mini Vi- 

As a consequence of the Corollary^, the right eigenvector v, 
which also represents the invariant probability of the Markov 
chain associated with the stochastic matrix A, plays a role in 
the controllability of the system represented by A. Vector v 
also measures the (right) eigenvector centrality of the network 
associated with A [2, Chapter 7.2]. It holds that, the bigger Vi 
is, the more relevant or central the node i is in the network, 
and therefore the fraction can be interpreted as a measure 
of heterogeneity in the node centralities. 

Exploiting this interpretation of v, we have that a network 
where all the nodes have similar centrality has a lower 
heterogeneity index and, according to Corollary 1, it 

will be more difficult to contror. On the other hand, easy 
to control networks need to have high heterogeneity in the 
nodes centrality. 

The evaluation of the heterogeneity index is particu- 
larly simple if A is reversible. The concept of reversibility for 
stochastic matrices is given in [25]. 

Definition 2 (Reversible matrix): An irreducible stochastic 
matrix A G satisfying A = A^ is called reversible. 

Remark 3: It is worth noticing that 


** Note that, with appropriate modifications, the Corollary can be applied 
also to row-stochastic matrices. 

^Interestingly, there is a connection between having a small heterogeneity 
index and the concept of wisdom of a network, as defined in [26], 


A := Cdmg{lZC)-\ (13) 

If C = then matrix A is reversible®. The latter matrix 
represents the transition matrix of the (weighted) random walk 
built on the network represented by the original matrix. 

Given a generic irreducible matrix C (possibly not symmet¬ 
ric), the right leading eigenvector of A, obtained from C as in 
(13), gives the (right) eigenvector centrality of the network 
associated with A, but it is also related to the PageRank 
[2, Chapter 7.4], which is an important centrality measure 
concerning in this case the network associated with O'. 

Assuming now that C is also symmetric, interesting analytic 
calculations can be carried out. Building A as in (13), its 
leading right eigenvector v has the following entries 

^col 

Vi = i = (14) 

with Cf°' := Cji and Ctot := J2i,j Consequently, the 
heterogeneity index is simply 

t^max ^ max,{Cf°^} 

^min min,{C'“i}' ^ 

S'mct is the degree of the node i in the undirected graph 
Q = {V,E,C), Formula (14) says that in this case the previ¬ 
ously defined centrality coincides, up to a normalization, with 
the so-called “degree centrality” of the network represented 
by C [2, Chapter 7.1]. 

Assume now that C is irreducible and symmetric and let 

c :=min{C'y | Cij 7 ^ 0}, 
c := max{Cij \ Cij 0 }. 

^Notice that, even if C is symmetric, A is not necessarily so. 

^Indeed for a generic matrix C the PageRank is defined as the right leading 
eigenvector of the column-stochastic matrix A = Q:Cdiag(lJC)“^ + (1 — 
a:)lnlX’ 0 < Q; < 1. The factor a., introduced to guarantee connectivity of 
the network, can be set to 1 since in our case C is irreducible by assumption. 
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Moreover, assume that all the diagonal elements of C are zero, 
i.e., Cii = 0, for all i. In this case, it turns out that 

W ^ {n-l)c ^ 

^min Q 

This shows that for reversible matrices obtained starting 
from matrices C for which c/c stays bounded, the fraction 
fmax/fmin Can gTow at most linearly in n. 

V. Examples 

In this section, we show how the previous reasonings 
can be applied to three examples. The first two examples 
concern random graph models, while the third one involves a 
structured graph. Since these examples will deal with matrices 
with zero diagonal, in order to avoid unnecessary technical 
complications*, as in [25] we use the concept of lazy version 
of a stochastic matrix. Precisely, from a column-stochastic 
matrix A, we can define a class of column-stochastic matrices 
as follows 

Aa '■= (1 ~ 0')A + Ctin , (17) 

where 0 < a < 1. Notice that for any a, A^ has the same 
leading right eigenvector of A and it is always primitive in 
case A is irreducible. 

Example 1 (Weighted random walk on a B-A graph): 
The Barabasi-Albert preferential attachment graph (B-A for 
short), denoted in what follows by BA(n, d), with d > 2, 
is a well-known and widely studied random graph model [2, 
Chapter 14]. Any realization of BA(n, d) is an undirected and 
unweighted graph with n nodes and, in its construction, each 
newly added node makes d connections with the previously 
existing nodes, according to a given rule that describes to 
which nodes the new node will be connected to (see [27, 
Section 4.1]). An example of a graph constructed using this 
model can be found in Figure 1. 



Fig. 1. realization of a B-A graph BA(50, 2) (the darker the node, the earlier 
it was inserted during the graph construction). 

In the following, we study the controllability degree of a 
lazy random walk on a weighted version of such a model. 

*In particular to avoid the presence of an eigenvalue in —1 and to also 
ensures that the product AA^ is primitive. 


Specifically, given an adjacency matrix (7 of a B-A graph 
constructed as above, we associate with it a symmetric matrix 
C obtained from C by letting Cij = 0 if Cij = 0 and 
otherwise letting Cij and Cji equal to a number drawn 
uniformly and independently at random in the interval [a,b], 
0 < a < & < oo. Notice that the matrix C defined in such a 
way is irreducible since the graph is connected. Using C, we 
can obtain matrix A as in (13), and fixing a constant a €]0,1[ 
we can build matrix Aa as in (17). The latter represents a lazy 
weighted random walk on a B-A graph. In order to obtain 
information on the controllability degree of the network using 
Theorem 1, we need to study the asymptotic behavior of a 2 - 
Using Cheeger’s inequality, it can be shown that the second 
eigenvalue of A, denoted with A 2 (A), satisfies 

A2(A) < B 

with high probability (w.h.p.) as n —> oo (see Appendix B for 
the details), where B ejO, 1[ is a constant depending only on 
a and b. On the other hand, the n-th eigenvalue A„(A) of A 
satisfies 

A„(A) > -1. 

These two facts imply, by (12), that® 

CT2 = max{Xl{Aa),Xl{Aa)} 

< max {((1 — a)z-f a)^} 
ze[-i,B] 

= max{((l — a)B + a)^, (2 q; — 1)^} < 1. (18) 

Thus, w.h.p. as n —> oo, the number a 2 is bounded by a 
constant smaller than one, depending only on a, a and b. Now, 
since by (16) Umax/^min < (ti — ^)b/a, Corollary 1 yields 

A{Aa,m) < -r. 

a <72(1 — (72 j 

By taking logarithm of both sides in the latter expression and 
making further computations we can argue that 

Tl 

log A{Aa, m) < k logn -log(l/(72), (19) 

m 

with A: > 0 being a constant depending only on a, b and (72 ■ 
Hence, since CT 2 is bounded away from 1 w.h.p. as n —t oo, 
the previous inequality allows us to conclude that, if m = 
o(n/log n), then 

logA(Ac,TO) _(x) w.h.p. 

or equivalently 

A{Aa,m) 0 w.h.p. 

This implies that the weighted lazy random walks on almost 
all the realizations of B-A graph model are difficult to control 
by means of o{n/ log n) driver nodes, as the cardinality of the 
network tends to infinity. 

Example 2 (Weighted random walk on a E-R graph): 
The Erdos-Renyi (E-R for short) graph model, denoted by 
ER(n,p), is one of the most celebrated and fundamental 
random graph models [2, Chapter 11]. A realization of the 

*Note that if A is an eigenvalue of A, (1 — a)A -f ct is an eigenvalue of 
Aa- 
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E-R graph model is an undirected and unweighted graph 
constructed as follows: starting from a graph of n nodes, 
we place an undirected and unweighted edge between each 
distinct pair of nodes independently and with equal probability 
p. It can be shown that [27, Theorem 2.8.1], if p = clogn/n, 
c > 1, the realizations of a E-R graph model are connected 
w.h.p. as n —> oo. With a partial abuse of language, we will 
refer to this subclass of E-R graph model as connected E-R 
graphs. Eigure 2 shows an example of a connected E-R graph. 



Fig. 2. A 100 nodes connected E-R graph. 

Using the same procedure described in the previous exam¬ 
ple, we build a symmetric matrix C with non-zero entries in 
[a, h] from the adjacency matrix C of a connected E-R graph. 
Erom the analysis given in [27, Section 6.5] and following the 
same lines of Example 1, we can argue that for the weighted 
lazy random walk on the graph described by C, the number CT 2 
is bounded away from 1 w.h.p. as n tends to infinity. Hence, 
we are in position to conclude that, by virtue of Corollary 1, 
the weighted lazy random walks on almost all connected E-R 
graphs are difficult to control w.h.p. as n goes to infinity if 
we use m = o(n/log n) control nodes. 

Remark 4: Interestingly, the choice made for the probability 
p allows us to make, in this case, further considerations on 
the heterogeneity index Wmax/r'min- As a matter of fact, for 
an unweighted E-R graph with p = clogn/n, it holds that 
the maximum and minimum degree (denoted by deg^^^ 
deg^in, respectively) satisfy [28, Ex. 3.4] 

degniax-^ 7 log n w.h.p. 

degmin-^ 7 log n w.h.p. 

where 0 < 7 < 7 < oo depend only on c. Consequently, since 
by (15) we have that 

/W ^ deg^ax b 
Vmin ~ degjai„a’ 

we can argue that Cmax/'^min is w.h.p. upper bounded by a 
constant as n —> 00 . Erom this fact we can conclude that a 
stronger result holds, namely that this class of systems are 
difficult to control w.h.p. as n goes to infinity, even if we use 
m = o(n) control nodes. 

Remark 5: The symmetry of the matrix C in Eormula (13) 
is necessary to enable a mathematical proof of the results 
proposed in the previous examples. In fact, if C is not sym¬ 
metric, we have no mathematical instruments to estimate CT 2 


and Umax/fmin- One could wonder how much this condition is 
crucial to observe practical uncontrollability of these networks. 
To investigate this issue, we carried out some simulations in 
which we did not impose symmetry on C. Eigure 3 shows that 
for B-A random graphs (T 2 stays bounded away from 1 and 
t^max/^min grows at most linearly in n, while for E-R random 
graphs CT 2 stays bounded away from 1 and Wmax/t-’min stays 
bounded. Erom this one can infer that symmetry seems not to 
be a crucial property and that the systems built in this way 
remains difficult to be controlled. 

Example 3 (Weighted random walk on a cube): Consider the 
unweighted graph consisting in the 3-dimensional fc-ary array 
[30], defined as the Cartesian product of three path graphs of 
length k (see Eig. 4). In analogy with the previous examples, 
starting from the adjacency matrix C of this graph, we build 
a symmetric matrix C having non-zero entries in [a, &]. Let 
A G M"^", n := k^, be obtained from C as in the previous 
examples and denote by the lazy version of A. 



Fig. 4. 3-dimensional 4-ary array: corner nodes are filled in black, while 
internal nodes in gray. 

By (15), we have that 

t^max ^ 26 

Wmin “ a 

since the maximum degree for the unweighted graph is 6 (for 
internal nodes) and the minimum is 3 (for corner nodes). In 
order to estimate A 2 (A) we exploit again Cheeger’s inequality, 
which yields that (see Appendix B) 

A 2 (A) < 1 - =: B{n), (20) 

where AT is a positive constant depending only on a and h. 
Now, using (18), it follows that 

CT2 = inax{Xl{Aa),XliAa)} 

< max{((l — a)B{n) + a)^, {2a — 1)^}. 

Notice that, for n sufficiently large, it holds that 

max{((l — a)B{n) + aY, (2q; — 1)^} = ((1 — a)B{n) + a^. 

The latter inequality can be used in (8), obtaining (see Ap¬ 
pendix B for detailed calculations) 

log(A(AQ,TO)) < ko- ki -h A:2logn 

m 
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Fig. 3. The asymptotic behaviour of <72 and fmax/nmin for networks with asymmetric matrix C, a = Xj'l and a = 1/2, b = 4. We start from an undirected 
random graph and, we associate with each edge (i,j) of the graph a weight Cij drawn uniformly and independently at random in the interval [a,b]. The 
upper figures refer to E-R graphs ER(n, 4 logn/n), while lower figures to B-A graph BA(n, 2). For each n, values are obtained averaging 250 realizations 
of each random model. The simulations were performed using Python’s NetworkX library [29]. 


where ko, fci and ^2 are positive constant values which depend 
only on a, b and a. When n tends to infinity we have 
that, as long as we use a number m = o(n^/^/logn) of 
control nodes, the system is practically uncontrollable since 
log(A(2lQ, to)) —>■ — oo, and therefore A{Aa,m) —)■ 0. 

Remark 6: In a similar way as Example 3, it is possible 
to prove that, given d > 3, the d-dimensional k-ary array is 

practically uncontrollable with a portion m = o ^ ”og ^ ^ of 

controllers. In case d = 1 or d = 2, our bound is not useful 
to determine the controllability degree of the system. 

VI. Conclusions and future directions 

In this paper we have shown the relevance of the network 
centrality for understanding how much energy is needed to 
control a dynamical network with adjacency matrix which is 
irreducible and (marginally) stable. Precisely, it is shown that 
if the matrix describing such a linear dynamical network has 
rather uniform node centralities, then it will be difficult to 
control. 

There are still many open questions that need to be ad¬ 
dressed. One is related to whether or not also the converse 
of the result shown in this paper holds true, namely if a 
system characterized by nodes with very different centralities 
is easy to be controlled. Simulations seem to disprove this 
conjecture, since node centrality alone seems not to be enough 
for characterizing the network controllability. 

Moreover, another interesting aspect concerns if node cen¬ 
tralities could be helpful in selecting control nodes, in case we 
want to increase the controllability degree. Indeed, one could 
argue that it is more convenient to control the nodes with 
largest centralities. Figure 5 contains some results obtained 
carrying out some simulations in order to clarify this question. 
In these figures it is shown the value of Aniin(VV’K:,T) as 
a function of n in three cases: random positioning of the 


control nodes, control nodes positioned at the nodes with 
largest centrality and control nodes positioned at the nodes 
with smallest centrality. Both for E-R and B-A random graphs 
something quite unexpected happen. Indeed in both cases 
positioning the control nodes at the nodes with smallest 
centrality happens to be the best strategy. We are unable to 
interpret this counterintuitive behavior'®, whose understanding 
surely deserves our attention in our future work on this subject. 

Another still open question is to prove what it seems quite 
intuitive, namely that if we can control a fixed percentage 
of the nodes, the energy needed to control the system stays 
bounded, or, more formally, if to = an, where a is any number 
in ]0,1[, then A{A,m) stays bounded away from zero. 

These are only a few questions that need to be better 
understood on this subject which, although it has received a 
lot of attention in the last years, still remains very challenging 
and full of issues that require to be investigated. 

Appendix A 

Proofs and additional results 
Proof of Lemma 1: We first note that the positive vector 

r := = [y/viWi ••• y/vfwfY 

is such that A^r = A^r, as can be seen by direct computation. 
Since A^ is a non-negative matrix and r is positive, we can 
use [21, Corollary 8.1.30] to conclude that A^ is the spectral 
radius of A^. Therefore ai = Xf ■ 

Proof of Proposition 1: We first prove two instrumental 
lemmas. 

Lemma 2: If x G M" is such that x^w = 0, then 

< '^2||a:||n, 

where CT 2 is the second largest eigenvalue of 

^^Note that there seems to be a connection with the results found in [4] 
concerning which are the most important nodes to control. 
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Fig. 5. Plots of the controllability degree Aniin(Vy;C,n) for B-A and E-R networks with symmetric matrix C, a = 1/2, a = 1/2, 6 = 2 and using m = n/3 
controllers selected in three different ways. Precisely, the controllers are placed at the nodes with highest PageRank centrality (Highest Centrality Nodes = 
HCN), with lowest centrality (Lowest Centrality Nodes = LCN), and at nodes selected at random (Random Nodes = RN). For each n, values are obtained 
averaging 500 realizations of each random model. 


Proof: First note that 

WA^xWfi = x^{A^)^nA^x = x^n^^'^A^U^/'^x. 

Since x^w = 0, we also have that is orthogonal to r, 

which is the eigenvector of A^' with respect to cri. Therefore 

x^n^/^A^n^/^x < a^x^VAl'^A^I'^x = (J 2 \\x\\^ 

and this ends the proof. ■ 

Lemma 3: If y G M" is such that y^v = 0, then 

Proof: Defining x := fL~^y, we have these two facts: 

x^w = y^fl~^w = y^v = 0 , 
l|a;||n = a;^na; = y^n-^mf-^y = ||2/||n-i- 
Moreover 

\\A^x\\l = x^nAn-^un-^A^iix 
= y^AU-^A^y 
= \\A^y\\u-^- 

Exploiting the previous lemma we have 

P^y|ln -1 = WA’^xW^ < <J2\\x\\1y = CTallylln-i- 

and we are done. ■ 

Using the lemmas we have just proven, we can finally prove 
the proposition as follows. Since y"' A*'~^v = = 0, 

defining x := {A^y-^y, we have that x^v = 0. Then 

\\{A^fy\\l.^=y^A^n-\A^Yy 
= x^AfJ-^A^x 

= P^^lln-i 

< a2x"'f\~^x 

= a2\\{A-^Y-^y\\l-. 

and applying the same reasoning t times, we obtain 

\\{A^Yy\\l..<4\\y\\u-^- 


Appendix B 

Detailed Calculations for Examples 1 and 3 


Eirst observe that, due to the Cheeger inequality [27, Section 

6 . 2 ] 


^ 2 (^ 1 ) < 1 


2 ’ 


( 21 ) 


where /i > 0 is known as the bottleneck ratio or Cheeger 
constant of the column-stochastic matrix A. This constant is 


defined as 


h = mm 
Scv 

v{S)<l/2 


Q{S,S) 

viS) 


( 22 ) 


where S' is a subset of V, S is its complement, v{S) := 
Efces^fc’ and 


Q{S,S):= Y. 

ies,jes 

Since in our case the column-stochastic matrix A is defined 
from a symmetric matrix C as in (13), then 

Q(S,S) = ^ ^ C,, >^e(S,S) (23) 


where e(S, S) denotes the number of edges between S and S 
in the unweighted graph described by C. The denominator of 
(22) can be treated as follows 


v(S) = 

kGS 


/^col 

kes 


L'tot 


(24) 


where vol(S) is the so-called volume of S in the unweighted 
graph and coincides with the total number of edges of the 
nodes belonging to S. Erom the previous inequalities we can 
argue that 

a . e(S,S) 
h> - mm —. 
b Scv vol(S) 

v{S)<l/2 


In [27, Section 6.4] it is shown that, for the B-A random 
graph, the minimum in the previous formula is lower bounded 
by a positive constant w.h.p. as n —> 00 . This shows that also 
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h is lower bounded by a positive constant w.h.p. as n —> oo. 
This solves Example 1. 

As far as Example 3 is concerned, we still use the Cheeger 
bound (21) and the inequalities (23) and (24). Then, since in 
this case vol(S') < 6 [S'!, we can argue that 


Q{S, S) ^ a e{S, S) 
v{S) - 6b 151 


which implies that 


, a 
h> — 
- 6b 


mm 

Scv 

\S\<n/2 


e{S,S) 

1^1 


where the minimum of e(5, 5)/|5| is the so-called isoperi- 
metric number of the graph [31, Chapter 3], which in a 3- 
dimensional fc-ary array is lower bounded by [30]. 

Therefore in this case we have 

a 1 
~ 3b v}!'^ 

Now, by virtue of the Cheeger inequality (21), it holds that 


A 2 (A) < l-iTn-2/3 =: B(n), 

2 

where K := is positive. 

Concerning the evaluation of A{Aa, m), by Corollary 1 and 
inequality (8), it follows that 


A{Aa,m) < 


2b B{n)^ 

~a B{n)^ (1 - B{n)) ’ 


where 


B{n) := ((1 — a)B(n) + a)^ = (1 — (1 — a)Kn 


Now take logarithm of both sides 

2& Ti — 

log(A(Aa, m)) < log-^-log B{n) 

a m 

— 2\og B{n) — log(l — B{n)). 


Using the Taylor’s expansion as n —> oo we get 

log i3(n) = 2 log(l — (1 — a)Kn~^^^) 

~ -2(1 -a)Xn-2/3 


and 

log(l - B{n)) = log(2(l - a)Kn-‘^^^ - (1 - 
~ log(2(l — a)Kn~‘^^^) 

= log(2(l - a)K) - ^ logn. 

Considering finally that —21ogi3(n) < 4(1 — a)K we can 
argue that 

\og(A(Aa,m)) < fco - ki -h fc 2 logn, 

m 

where ko, fci and fc 2 are positive constant values which depend 
only on a, b and a. 
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